Synchronization Defects and Broken Symmetry in Spiral Waves 
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Spiral waves are investigated in oscillatory media exhibiting period-doubling bifurcations. In the 
period-doubled and chaotic regimes, the rotational symmetry of the spiral wave is broken. The loss 
of symmetry takes the form of synchronization defect lines where the phase of the local oscillation 
changes by multiples of 2tt. The internal structure and geometry of these synchronization defects is 
studied and a discussion of the possible types of defect lines is presented. 



ON 

o 

Q 



> 

o : 
o 

ON . 

"a : 

o : 

^ . 
o ■ 



■ 

X 



82.20.Wt, 05.40.4-j, 05.60.+W, 51.10.-fy 



Spatially-distributed oscillatory media may undergo 
bifurcations where the period of the orbit doubles at al- 
most every point in the system. In the simple oscillatory 
regime (period-1), before such period-doubling bifurca- 
tions occur, the system may support stable rotationally- 
symmetric spiral wave solutions, jl] The dynamics in this 
regime is typically described by the complex Ginzburg- 
Landau (CGL) equation, the generic equation for an 
oscillatory medium near the Hopf bifurcation point at 
which oscillations appear Spiral waves are a well- 
known feature of this regime where they have been in- 
tensively studied However, they can also exist in 
media where the local dynamics is complex-periodic or 
chaotic [Q-^. In this Letter, we investigate the conse- 
quences of period-doubling bifurcations on the structure 
and dynamics of spiral waves. We find that the symme- 
try of the spiral wave is broken by defect lines where the 
phase of the oscillation changes by multiples of 27r, and 
we study the nature of these synchronization defect lines. 

We consider the dynamics of spatially-distributed sys- 
tems governed by reaction-diffusion equations of the form 



dc{r,t) 
dt 



R(c(r, t))+DV'^c{v, t) 



(1) 



where c(r, t) is a vector of local concentrations, D is the 
diffusion coefficient and R(c(r,t) describes the local re- 
action kinetics. While the phenomena we investigate 
should be observable for any spatially-distributed sys- 
tem exhibiting period-doubling, we have considered cases 
where the spatially- homogeneous system c{t) = R(c(t)) 
itself exhibits period-doubling bifurcations. Specifically, 
the calculations described here were carried out on the 



Rossler model, Rx = —Cy — c^, Ry = 



ACy , Rz 



CxCz — Ccz + B, which is well known to exhibit chaos 
arising from a cascade of subharmonic bifurcations. 

Spiral waves were initiated as in [^,|| , taking advantage 
of the cyclic character of the projection of the Rossler at- 
tractor on the {cx, Cy) plane. Various values of the param- 
eter C in the interval [2.5,6.0] were considered, while the 
other parameters were fixed a.t A = 0.2 and B = 0.2. The 
scaled diffusion coefficient was DAt/ (Ar)^ = 1.6 x 10"^ 
{At = 10"^) in all calculations. Simulations were car- 
ried out on a disk-shaped domain of radius 255 with no- 
flux boundary conditions. In the period-1 regime, the 



medium supports a single, stable, one-armed spiral wave. 
As the parameter C increases, the local dynamics under- 
goes a period-doubhng bifurcation at C = C* ~ 3.03. 
(The spatially-homogenous Rossler model bifurcates at 
^odc ~ 2.83). In order to quantitatively characterize 
this transition, it is convenient to introduce a local order 
parameter i5c2(r), defined as the difference between two 
successive maxima of the time series C2(r,t). 
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FIG. 1. Bifurcation diagram, constructed at ro = 130, 
showing the first period-doubling bifurcation in the medium 
supporting a single spiral wave (diamonds). For comparison, 
the upper left curve is the bifurcation diagram for a spatially 
homogeneous system. 

Figure]^ shows the dependence of Scz{ro) on C calcu- 
lated at a typical point at radius rg — 130, sufficiently 
far from both the boundary and the core of the spi- 
ral. One observes a parabolic increase in Scz{ro) beyond 
C* ~ 3.03, typical of a supercritical period-doubling bi- 
furcation. This behavior is seen at almost every point in 
the medium (the structure is difhcult to resolve deep in 
the core region). Along almost any ray emanating from 
the core, the magnitude of Scz{r) varies as A{r)^/C — C* 
where A{r) behaves like the amplitude profile of the spi- 
ral itself, i.e., A(r) ~ in the core region and is constant 
elsewhere except near the boundary. 

Within the period-2 domain, the spiral wave acquires 
a global structure different from that in a simple peri- 
odic medium. Figure ^ shows the Cz{Y,tQ) concentration 
field at a single time instant t^. The alternation of high 
and low Cz(r,to) maxima unambiguously demonstrates 
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that period-2 local temporal dynamics manifests itself in 
the formation of a period-doubled spiral waveform with 
broken rotational symmetry whose wavelength is twice 
that of the original spiral wave in the period- 1 medium. 
The lower panel of Fig. || shows Cz(r,to) in grey shades 
and indicates the curve connecting the spiral core and 
the boundary, denoted as Q, where sharp changes in the 
concentration occur. The curve plays a central role in 
the organization of the spiral wave, and its character will 
be examined below. 




FIG. 2. Spiral wave in a medium with period-2 local dy- 
namics at C = 3.84. Concentration field Cz{r,to) is shown as 
elevation in the upper panel and as grey shades in the lower 
panel. The solid line depicts the Q curve. The arc segment at 
radius rn ~ 130 along which points were taken to construct 
Figs. |^-|5| is also shown. 

The fact that the local dynamics is almost everywhere 
period-2 endows the spiral with some unusual features: 
unlike a period- 1 medium where the concentration is pe- 
riodic with the period of the spiral rotation, here, after 
one turn of the spiral, the high and low maxima inter- 
change and it is only after two spiral periods that the 
concentration field is restored to its initial value. Al- 
though the spiral rotates, the curve is, up to numeri- 
cal accuracy, stationary. In the asymptotic regime, after 
transients depending on the initial condition, the shape 
of the curve takes the form of a straight line segment 



with a short curved portion lying inside the core region. 
Figure § shows 6cz {0) for three values of C as function of 
the polar angle 9 along the arc segment that crosses the 
curve. To a very good approximation, the numerical 
data show that Scz{r, 9, C) varies like 



Scz{r,9,C) ~ A{r)VC-C* taTih[Kr{9-9n)iC-C*)] (2) 

where 9q is the angular position of the fl curve and k is a 
numerical factor. The Q curve is thus a localized object 
whose width remains approximately constant with r and 
varies like 1/(C — C*). As C approaches C* from above, 
the width increases until it becomes comparable to 27r 
and the H. curve ceases to be a well-defined object [Q. 
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FIG. 3. Order parameter &Cz[Q) profile across the 57 curve 
along the arc at radius ro = 130 indicated in the lower panel 
of Fig. I for three values of C: (a) C = 3.04, (b) C = 3.10, (c) 
C = 3.84. Curve (c) is calculated for the medium in Fig. y. 

We now show that the curve corresponds to a 
one-dimensional synchronization defect across which the 
phase of oscillation changes by some multiple of 27r. One 
can always uniquely parametrize an orbit of a period- 1 os- 
cillation by a phase variable ip S [0, 27r) and describe the 
spatially-distributed medium by a phase field (p(r, t). At 
the center of the spiral lies a point defect of this field char- 
acterized by a topological charge ^ § V(^(r, i) • dl = nt 
1^, where the integral is taken along a closed curve en- 
circling the defect. There exist several possible ways to 
extend the definition of phase for complex-periodic os- 
cillations 1^. We assume that the phase of an n- 
periodic oscillation is a scalar function of time which in- 
creases monotonically by 27rn for each full period of the 
oscillation T„. For some systems, can be defined in 
terms of an angle variable in a suitably-chosen coordi- 
nate frame in phase space. For the Rossler model we 
use a cylindrical coordinate system (p, c/?, z) with origin 
at the unstable fixed point of the spatially-homogenous 
system and z along c^- The phase variable G [0, 2?™) 
then takes the form = + 27rm where m is an integer 
with values from to n — 1. While is a single- valued 
function of the original dynamical variables, the phase 
of a complex-periodic oscillation is not an observable, 
since a knowledge of the entire local orbit is required in 
order to calculate to. 
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FIG. 4. Two Cz(i) concentration time series calcu- 
lated for the medium shown in Fig. ^ at points 
ri,3 = (ro — 130, 6*^ ± 0.05) on opposite sides of the Q curve. 

The n curve can be viewed as a one-dimensional defect 
of the (/)(r, i) field. Figure ^ shows two Cz(t) concentra- 
tion time series at nearby spatial points on either side of 
the O curve. The two oscillations are shifted relative to 
each other by half a period. For a period-2 oscillation, 
this corresponds to a phase shift of Scj) = I-k. The neces- 
sity for such a phase synchronization slip in a one-armed 
spiral in a period-doubled medium can be understood 
from the following argument. Consider a contour inte- 
gral of the phase gradient V(/)(r, t) taken along a closed 
loop r which surrounds the core of the spiral wave. For 
an arbitrary point rg G F, this loop is just a path in 
the medium which starts at Tq, described by the local n- 
periodic oscillation phase (/)(r, i), and returns to the same 
point with same oscillation phase (f)(r,t). Thus, the in- 
tegral may take only values 2'Knk equal to multiples of 
the full-period phase increment. Since the phase field is 
given by (/>(r, t) = ip{r, t) + 27rm(r, t), we find 

^V0(r,i)-dl= <j>\7ip{r,t)-d\ + 2'K <j>\7m{r,t)-d\ . (3) 

Simulations for a one-armed spiral show that, at any 
given time t, § V(/9(r, t) ■ d\ — ±2tt, regardless of the pe- 
riodicity of the local dynamics |l^. Given that the full- 
period phase increment is 2tt only in the case of period- 1 
dynamics with m{r, t) = 0, for period-doubled dynamics 
the integration of Vm(r, i) along F must yield a non-zero 
contribution to balance (j^). Since m(r,t) is an integer 
function, its value changes discontinuously with time and 
space so that Vm(r, t) is different from zero only at a sin- 
gle point: the intersection of F with the O curve. For a 
period-2 medium, the addition of a 2tt phase jump on the 

curve and a 2tt contribution from integration of ip{r, t) 
yields the necessary full-period increment of Att. 

The nature of the phase jump associated with the f2 
curve can be understood from the observation of the lo- 
cal orbit loop exchange which occurs on the scale of the 
il curve width. Figure H shows local orbits calculated at 
ro = 130 for three angles Oi = Oq - 0.05,6*2 = On, 83 = 
6n + 0.05. (The corresponding Cz{t) time series for the 
first and last orbits are shown in Fig. ^ and 6cz [0) is given 
in Fig. ^(c).) Consider local orbits on the arc connect- 
ing points 61 and ^3 (cf. Fig. ||). As one traverses the 



arc, the larger, outer loop of the local orbit constantly 
shrinks while the smaller, inner loop grows. At ^ = ^n, 
both loops merge and then pass each other exchanging 
their positions in phase space. (Compare Figs, ^(a) and 
(c).) The behavior of the 5cz{0) order parameter along 
the arc provides only limited information on the loop ex- 
change. As one sees from Fig. ^ not only are the Cz 
maximum values of two loops equal {5cz{9n) = 0), but 
the entire loops coincide in phase space: at the exchange 
point 9 = 9q, the local oscillation is effectively period- 1 
. Continuity of the medium requires that the period- 1 
points form a line extending from the core to the bound- 
ary: the Q curve. 




FIG. 5. Loop exchange in local orbits as the curve is 
crossed (see text). Data collected on the arc shown in Fig. ^: 
(a) 01 = 9n- 0.05; (b) 62 = 9n; (c) 83 = 9n + 0.05. Two 
points corresponding to time instants ti = 10.3 (diamond) 
and t2 = 19.8 (asterisk) are marked on each trajectory to 
highlight the exchange. 

The n curve separates domains that are dynamically 
equivalent apart from a phase shift and it should have no 
net velocity in its normal direction. Outside the core re- 
gion, where the Q curve width is small on the scale of the 
spiral wavelength, any large-scale curvature will be elim- 
inated by motion of the f2 curve. This motion, induced 
by the mean curvature and proportional to the diffusion 
coefficient, will yield a straight curve. In the core re- 
gion, the curve width is comparable to the length scale 
of other concentration gradients and more complicated 
structure is possible. 

We now consider more complex local oscillations. At 
C « 4.09, the system undergoes a bifurcation to period 
4 which can be described by an order parameter 52Cz(r), 
analogous to 6cz{r), defined as the difference between 
the first and third or second and fourth Cz(r,t) max- 
ima within one full period of the oscillation. This can 
be generalized to any period-2'^ regime. An oscillatory 
medium with period-4 local dynamics has two types of 
synchronization defect lines: Oi and Q2 curves with as- 
sociated phase jump values of ±67r and An, respectively 
|]lT| . These curves are shown in Fig. ^ for C — 4.3. 

We may number the loops of the local period-4 cycles 
successively according to their positions in phase space 
starting from the innermost loop. Then the fli curve 
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corresponds to the (43^!) loop exchange altering all four 
loops, while the curve is attributed to the (2143) ex- 
change which involves only rearrangement of loops in- 
side the period-2 bands. By considering topological braid 
properties of a period-doubled orbit [p| , all possible loop 
exchanges can be enumerated. Outside the core region, 
we expect that there exist k types of curves in a period- 
2^^ medium. 




FIG. 6. Synchronization domains delimited by Sli and 0,2 
curves in a medium with period-4 dynamics. The maximum 
value of Cz(r) in one spiral rotation is plotted as grey shades. 
Apart from boundary effects, roughly four levels are observed, 
corresponding to the four loops of the local orbit, from loop f 
(darkest shade) to loop 4 (lightest shade). 

The argument given above for the value of the 
/ V(/)(r, t) ■ dl integral applies to a period-4 medium as 
well. The existence of two types of defect provides sev- 
eral ways to satisfy condition (^. Thus, a medium with 
a one-armed spiral wave characterized by nt = +1 may 
have one fli curve with -t-67r phase shift or both fli (— Btt) 
and 0,2 (47r) curves. Configurations where more than two 
curves originate in the same spiral core are, in principle, 
possible. However, they were not observed starting from 
the initial conditions considered here. 

For higher values of C, the Q curves may evolve in 
time leading to spatio-temporal chaotic regimes within 
the general spiral structure. The curves separate do- 
mains of near synchronization where the local phase 
(j){v, t) changes continuously but experiences sudden 
jumps when a domain boundary is crossed. A full ex- 
ploration of the dynamics of O curves in this and other 
parameter regimes remains to be carried out. 

Wavelength-doubled spiral waves have been observed 
in recent experiments on the Belousov-Zhabotinsky (BZ) 
reaction in ruthenium-complex monolayers on the surface 
of a BZ solution. |]l2| The experimentally observed spiral 
wave structure with a clearly visible Q curve (cf. Fig. 5 
of Ref. jl2| ) is identical to that in Fig. |2[ Not only is the 
final wavelength-doubled spiral wave sturcturc the same, 



but the dynamics leading to its formation from an initial 
period- 1 spiral is also the same as that observed in the 
experiment. We have demonstrated that such broken- 
symmetry spiral waves are necessarily formed when the 
local dynamics has period-2 character, as is the case in 
Ref. . Although our results were obtained for a system 
possessing a period-doubling cascade, we believe that the 
synchronization defect lines are generic features of dis- 
tributed media exhibiting complex periodic behavior, in- 
dependent of the specific origin of this periodicity. Thus, 
we expect the phenomena described in this Letter to be 
observable in excitable as well as oscillatory media where 
such period-n local dynamics is observed. 
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